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SUMMARY 


A physical interpretation of the observed form of the pressure distribu- 
tion beneath a two-dimensional "short" separation bubble (which modifies the 
external inviscid pressure distribution only locally) is given in terms of 
boundary layer concepts (i.e., constancy of static pressure across the layer as 
long as the layer is "thin"). At the mean separation and reattachment points 
(which lie on the same mean streamline) , the local static pressure equals the 
local stagnation pressure, since the velocity is zero at these points. The 
boundary layer hypothesis then implies that reattachment can only occur at a 
point X - Xg downstream of the separation point Xg if the jump in external 


inviscid static pressure 


-[u2(x) - u2(Xs)J 


is less than or equal to the rise 


in stagnaticn pressure H(x) - H(Xg) along the separation streamline after 
separation. A simple method for estimating the growth of H(x) along a mean 
streamline entrained into the underside of a growing shear layer through the 
transition region is discussed, and predictions of bubble bursting conditions 
and a lower bound on the bubble length are compared with experiment. 


INTRODUCTION 

Since the early work of Melvill Jones (ref. 1), it has been known experi- 
mentally that a laminar boundary layer on an airfoil, after entry into a region 
of adverse pressure gradient strong enough to cause laminar separation, will 
separate in the laminar state, sometimes achieve transition to turbulent flow 
in the separated shear layer, and (if the adverse pressure gradient is not too 
severe) reattach to the surface to form a closed recirculating flow eddy known 
as a separation bubble. In a typical case (fig. 1), the turbulent boundary 
layer downstream of the reattachment point R either does not separate at all 
over the remaining portion of the upper surface or separates only a short dis- 
tance upstream of the trailing edge. In such a flow, the form of the external 
static pressure distribution outside of the region of separation is approximated 
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reasonably well by the potential flow pressure distribution, though the magni- 
tude of the actual pressure is somewhat smaller. It follows that the lift of 
an airfoil with this type of bubble separation is approximately equal to the 
value predicted by inviscid theory. Following Tani (ref. 2), we will refer to 
this flow phenomenon as a separation bubble of the "short" type, to distinguish 
it from another type of closed recirculating flow in which the entire pressure 
distribution differs radically from the potential flew form. The latter is 
referred to as a "long" bubble. 

Separation bubbles are of great importance in engineering because of the 
role they play in the phenomenon of airfoil stalling, which may sometimes be 
identified with a transition of the flow from a short type of bubble separation 
to the long type. This abrupt transition, known as "bursting" of the short bub- 
ble, is in some cases responsible for producing the well-known "critical" Reyn- 
olds number region in the plot of drag coefficient versus Reynolds number for 
circular cylinders and spheres (ref. 2) and for moderately thick airfoils at 
Reynolds numbers lower than about 10^ (refs. 3 and 4) . 

Reference 2 contains a thorough background discussion of the phenomenon of 
separation bubbles and a review of the experimental work prior to 1964. Refer- 
ence 5 reviews most of the more recent experimental work up to 1976. The reader 
is referred to these sources for general information about the subjects treated 
in this paper. 

In recent years, there have been several attempts to develop semi-empirical 
calculation methods for predicting the length and bursting conditions for sepa- 
ration bubbles of the short type. These range in sophistication from solution 
of the full time-dependent two-dimensional Navier-Stokes equations (with model- 
ing of the turbulence), as in reference 7, to simpler modifications of existing 
boundary layer calculation methods incorporating iterative methods to account 
for viscous-inviscid interaction between the external stream and the flow within 
the bubble (refs. 6 and 8) . All of these methods use experimental data on sepa- 
ration bubbles to fix the values of certain numerical coefficients appearing 
in the theory. In each case, the shape of the separation streamline and the 
external pressure distribution between separation and reattachment are unknowns 
in the problem and are not determined until the calculation is finished. Since 
the distributions of these quantities play a crucial role in determining the 
length and bursting of the bubble, there is some difficulty in using the numerical 
prediction methods (or the theories on which they are based) to deepen one's 
understanding of the fundamental physical processes controlling the flow. 

The objectives of the present study are to explain why reattachment occurs 
at all, to explain why the observed pressure distribution has its characteristic 
form, and to develop a shortcut method for calculating the length of a short 
separation bubble, if one exists, or to determine the conditions under which a 
short bubble cannot exist. 
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PHYSICAL MODEL 


The Cause of Reattachment 

The principal cause of turbulent reattachment following laminar separation 
has been discussed by Cebici and Bradshaw (ref. 9) . Consider the separation of 
a laminar boundary layer fron a smooth surface as shown in figure 2. The edge 
of the shear layer has been drawn to shew the spreading of the shear layer in 
the downstream direction due initially to laminar and later to turbulent mixing 
of momentum. This spreading results in a greater mass flux across the seg- 
ment BB' in figure 2 than across the segment AA‘ . In two-dimensional flow, 
mass conservation requires that there be a net inflow across the segment AB. 
This entrainment of fluid into the underside of a growing shear layer will be 
greatly increased by the transition from laminar to turbulent flow. If this 
entrainment rate is greater than any reversed flew coming from far downstream, 
then the shear layer simply sucks itself back onto the surface from which it 
separated (ref. 9, p. 366) to form the short separation bubble shown in fig- 
ure 3, in which the entrainment flow is supplied by the splitting of the shear 
layer at reattachment. 


The Form of the Pressure Distribution 

In trying to identify the most important factors controlling the dynamics 
of separation bubbles, it is useful to see how many of the qualitative features 
of the flow can be accounted for adequately by simple boundary layer concepts. 
If, for example, simple boundary layer theory is found to account adequately 
for the observed length and pressure distribution of a bubble, then one can 
infer that the direct effects of streamline curvature and nonzero pressure gra- 
dients in the cross-stream direction (which are neglected in boundary layer 
theory) are indeed unimportant, at least in determining the gross qualitative 
features of the flow. This process of testing hypotheses that certain terms in 
the equations of moticxi are negligible will serve to increase our understanding 
of what a separation bubble is if it allows us to identify the few important 
frown among the many unimportant effects operating at once in the flow. 

In this spirit, we attempt to explain the characteristic pressure distri- 
bution beneath a short separation bubble (cf. fig. 1) by supposing that the 
assumptions of classical boundary layer theory hold. In particular, the static 
pressure across the separation bubble in figure 3 is assumed to be constant, and 
the velocity component in the direction normal to the wall is assumed to be an 
order of magnitude smaller than the component parallel to the wall. With these 
assumptions, the pressure distribution along the wall between points S and R 
in figure 3 can be inferred. The pressure gradient between points S and K 
must be small since the velocities and shear stresses in the slowly moving 
interior of the bubble are driven only by entrainment and are, therefore, small. 
In the region downstream of the point K, the streamwise velocity along each 
closed recirculating flow streamline will be zero at the downstream extremity 
of that streamline. Since the vertical velocity at such a downstream extremity 
is small, the static pressure is very nearly equal to the local stagnation pres- 
sure at that point. It follows that the distribution of pressure along the line 
connecting all such downstream extremities of streamlines (represented as the 
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dotted line in fig. 3) will approximately equal the distribution of stagnation 
pressure between the lower edge of the shear layer (at point K) and the sepa- 
ration streamline (at point R) . This difference in stagnation pressures may 
be substantial because the velocity on the separation streamline (which lies 
close to the center of the shear layer at the streamwise station K) will typ- 
ically be much larger than the velocity on the lower edge of the shear layer 
at any given streamwise station. 

In general, we expect the pressure distribution beneath the separation bub- 
ble in figure 3 to be as shown in figure 4, which is in good qualitative agree- 
ment with the examples of short separation bubbles in reference 2. The rather 
sharp corner in the pressure distribution in figure 4 is interpreted here as 
the point where the lower edge of the shear layer impinges on the wall. Accord- 
ing to our model, this point also coincides with the center of the innermost 
closed streamline in figure 3. Note that the point K is not identified with 
the transition from laminar to turbulent flow (which we assume takes place some- 
where upstream of the point K) . 

It has been assumed in all experimental studies on separation bubbles known 
to the author that the sharp corner in the pressure distribution coincides with 
the transition point- Gault (ref. 10) proposes this as a means of experiment- 
ally determining transition, while Van Ingen (ref. 11) uses the same assumption 
to fix the numerical values of the adjustable constants in a form of his "e to 
the ninth" method for transition prediction, developed especially for use in 
separation bubbles. In view of the fact that transition from laminar to turbu- 
lent flew occurs quite readily in a free shear layer in zero pressure gradient 
without modifying the static pressure, one is at a loss to see why transition 
as such should have any effect on the pressure distribution at all, much less 
account for the abrupt, nearly discontinuous pressure rise that occurs at the 
rear of a short separation bubble. The rather arbitrary and ad hoc character 
of the conventional interpretation of the pressure distribution renders it less 
favorable from a fundamental viewpoint than the interpretation in terms of stag- 
nation pressures given above. 


Length and Bursting of Short Bubbles 

We have included among the defining features of short separation bubbles 
the condition that the bubble modify only locally the static pressure distribu- 
tion about the body on which the bubble occurs; that is, the pressure distribu- 
tion everywhere ahead of and behind the bubble is approximately equal to the 
value predicted by potential theory. In particular, the pressures at the sepa- 
ration and reattachment points (which are the inner limits of the portion of the 
pressure distribution lying outside the bubble) must both be points lying very 
near the inviscid pressure distribution curve. These endpoint conditions, 
together with our model of the pressure distribution within the bubble in terms 
of stagnaticxi pressures on the recirculating flew streamlines, provide a clue 
as to what determines the length of a short separation bubble, when it is pos- 
sible for one to exist, and why it is sometimes impossible for a short bubble 
to exist. 
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Now the difference in static pressure between the points S and R in 
figure 3 equals the net rise in stagnation pressure along the separation stream- 
line, since the velocity is zero at both endpoints. It follows that if a plot 
of the rise in stagnaticm pressure along the separation streamline downstream 
of separation is drawn on the same graph as the rise in inviscid static pressure 
downstream of the same point, then an estimate of the bubble length is the 
streamwise distance from separation to the first point where the curves cross. 

One difficulty with such a graphical procedure for the calculation of the 
bubble length is that the curves may cross at more than one point, giving a 
choice of several possible bubble lengths rather than a single one. Another 
difficulty is that the stagnation pressure along the whole length of the sepa- 
ration streamline is not known under general conditions. It can, however, be 
estimated by assuming that the growth of the shear layer downstream of separa- 
tion is nearly equal to that of a free shear layer in zero pressure gradient, 
for which the mean velocity and, hence, stagnation pressure and stream function 
profiles are known from experiments on transition in a laminar mixing layer 
(ref. 12). 

A calculation based on similar physical concepts was made by Tani (ref. 2) 
to estimate the maximum attainable value of the coefficient of pressure recovery 
within a bubble o defined by the relation 

Pr “ Ps 

a = (1) 



Owen and Klanfer (ref. 13) had proposed the value = 0.35 from experi- 

mental evidence. By noting that the velocity along any streamline entrained 
into the low velocity side of a growing shear layer cannot exceed the velocity 
on the dividing streamline (^/U)d.S. ~ 0.5873 ..., Tani argued that a lower 
bound on the nondimensional stagnation pressure rise existed and was given by 

J = [(u/U)d.sJ ^ = 0.3449 ... (2) 

max 



which is in excellent agreement with the figure 0.35 given by Owen and Klanfer. 
Tani did not suggest the extension of the stagnation pressure idea to the calcu- 
latiOTi of the bubble length, however. 

The method of estimating the reattachment position as the streamwise coord- 
inate of the point of intersection of the curves of inviscid pressure rise and 
stagnation pressure rise on the separation streamline after separation provides 
a necessary condition for the existence of a short bubble, namely, that the two 
curves indeed intersect. A possible condition for the impossibility of a short 
separation bubble (i.e., a "bursting" condition) might then be that the inviscid 
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pressure gradient downstream of separation be so steep that the two curves never 
intersect. The borderline case in which a short bubble is just barely possible, 
but would become impossible with an infinitesimal increase in the inviscid pres- 
sure gradient, is when the two curves are tangent to each other at a single 
point. Later in this paper we will use this tangency condition to derive a 
close approximation to Caster's "bursting line" (ref. 14), which has been found 
by several experimentalists (see ref. 5) to be a reliable empirical formula for 
predicting the boundary between the short and long type of bubble separations. 


DEVELOPMENT OF THE SHEAR LAYER 

For reasons indicated in the last section, it is desirable to have an 
explicit formula for the development of the stagnation pressure H(s) along 
the separation streamline in a short bubble. 

For simplicity, we will assume that the actual distribution of stagnation 
pressure along the separation streamline can be approximated by the stagnation 
pressure distribution along a particular streamline entrained into the underside 
of a self-similar shear layer in zero pressure gradient, for which the velocity 
profile has been calculated from the Blasius equation in the laminar case 
(ref. 15), and has been measured experimentally in the transitional and fully 
turbulent cases (refs. 12 and 16, respectively). It has been found experimen- 
tally (cf. ref. 12, fig. 12) that the functional form of the mean velocity pro- 
file is experimentally indistinguishable for these three cases, the effect of 
transition being a large and abrupt increase of the rate of spreading of the 
shear layer measured by, say, the growth of the momentum thickness 6 (s) . 

Before treating the problems of the location of the transition region, the 
growth of the resulting turbulent layer, and the determination of the value of 
the stream function (relative to the reference streamline in a self-similar 
shear layer) appropriate for representing the separation streamline in a sepa- 
ration bubble, it is necessary to recall certain properties of self-similar 
shear layers, and, in particular, to obtain explicit approximate formulas for 
the velocity and stream function profiles. 


The Laminar Shear Layer 

The nomenclature for self-similar shear layers is as shown in figure 5, 
which illustrates the flow downstream of a hypothetical splitter plate, the 
upper surface of which moves with the velocity U-| (precluding the formation 
of a boundary layer there) , while the lower surface is stationary. The fluid 
above the plate (x < 0, y > 0) moves with the uniform velocity U-| while the 
fluid below the plate (x < 0, y < 0) is at rest. The effect of these boundary 
conditions is the elimination of any length scale other than the streamwise 
coordinate x, so that self-similarity in laminar flow is to be expected. 

The steady two-dimensional Prandtl boundary layer equations for zero pres- 
sure gradient. 
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thus reduce to the Blasius equation, 


2 f ' " (n) + ff " (n) = 0 


in the usual way (e.g., ref. 17, p. 126), where 



vlJ(x,y) vKx,y) 

f(n) = = 

(VUix)V2 Ut il(x) 


34; 

34; 1 /Uiv\ 

"1 

u = — = f (n) 

3y 

3x ~ 2\ X / L 

n f (Ti) - f (n)J 


(4) 


(5) 


( 6 ) 


(7) 


( 8 ) 


The boundary conditions appropriate for a free shear layer are 


U ->■ 

Ui f ->■ 1 

as 

n ^ 

+ oo 

(9) 

U -> 

0^ f< -v 0 

as 

n ^ 

— oo 

no) 

4; = 

o 

II 

> 

II 

d.1 

when 

n = 

0 

(11) 


The exact solution of equation (5) satisfying the boundary conditions (9) , 
(10), and (11) has been tabulated by Lock (ref. 15) to four significant figures. 
The tabulated solution provides a useful check on the accuracy of approximate 
representations of the velocity profile in terms of elementary functions. 
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^^proximation of the Velocity Profile 

For the present purpose, it is sufficient to obtain an approximate solu- 
tion to equation (5) in terms of elenentary functions whose forms are chosen 
to reproduce as many of the properties of the exact solution as is convenient. 
For example, the condition f(0) = 0 implies that f * ' ' (0) = 0 (from equa- 
ticxi (5)), so that the velocity f has an inflection point at the origin, or, 
equivalently, the shear stress is maximum there. In addition, the conditions 
f(Ti) ~ n as n -»■ +00 and f(ri) « Constant as r| imply that f 

behaves like a complementary error function for large positive n and like a 
simple exponential for large negative p. These conditions, together with the 
requirement that the velocity profile be a smooth monotonic function of n, 
suggest the approximate forms 


n > 0; f ' (n) 


[’ 


£’(o^- r 

/it ‘^ari/2 



n < 0; 


f (n) = f (0) 


0 

- sech'^ 
2 


bp + tanh“^ 



( 12 ) 


(13) 


each of which has an inflection point at p = 0 and the correct asymptotic 
behavior as |p| Equations (12) and (13) are constructed so that the 

velocity f ' (p) is continuous at the juncture p = 0. The requirement that 
the slope of the velocity be continuous at p = 0 provides one equation relat- 
ing the constants a and b: 


-[1 - f • (0)] a 




(14) 


Equation (14) provides an expression for the ratio a/b in terms of f'(0), 
where either a or b remains free to be determined by a suitable 
normalization. 

The only important unknown that remains is now f'(0), which we will deter- 
mine by a momentum integral method. Integrating equation (5) frcxn -<» to 0, we 
obtain 


2 f 


+ ff 



(f')2 dp 


0 


while the integral fran 0 to «> gives 
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f (-1 + f) dn = 0 


2 f 



+ 


f (f ' 



Applying the boundary conditions for a shear layer , these equations beccane 


2 f • (0) 



(f )2 dn 



f {1 - f ) dn 


( 15 ) 


which are forms of the Von Karman momentum integral equation (e.g., ref. 17, 
p. 146) for each half of the flow, and which have direct counterparts when the 
flow is turbulent (e.g., ref. 18, p. 227). 

If equations (12), (13), and (14) are substituted into equation (15), then 
a cubic equation for the quantity 1 - f • (0) is obtained (see appendix) , the 
relevant solution of which gives f'(0) = 0.58923, which agrees well with the 
tabulated value f ' (0) = 0.5873 in reference 15. 

For our normalization condition (needed to determine either a or b from 
equation (14)), we note that equation (15) implies 


e 

- E 0 

z 


* 


C f (1 - f) dn 


^ oo 

+ \ f (1 - f) dn 
'^0 


1° 

= f - f" (0) + f " (0) = -f (-<») 


(16) 


which according to reference 16, has the numerical value -f (-<») = 1.2386. 


It follows from equations (12), (13), (14), and (16) (see appendix) that 
a = 0.88544 and b = -0.30160, which completely determines all of the constants 
in our approximate velocity profile. For n < 0 (which is the region of 
interest in the separation bubble problem), the approximate velocity profile 
given by equation (13) is uniformly accurate to within 1 percent of the tabu- 
lated values in reference 16. 


Development of the Momentum Thic)cness 

The nondimens ional velocity profile f ' (D) , together with the definition 
of the transverse length scale (x) in equation (6) and the free stream veloc- 
ity U-| completely determine the flow in a laminar shear layer. It is conve- 
nient to express Z (x) in terms of information available at a station x = x-| , 
rather than in terms of the distance x downstream of the splitter plate in 
figure 5. Prom equation (6) , 


185 



/vx\^/2 

V(x - xt) 

il(x) = — 

+ a (Xt) 2 

\Ul/ 

L ut j 


1/2 


= £■ 


1 X - XT 




+ 1 
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1/2 


where il-] is defined to be S,(xi). Rewriting in terms of the monentum thick- 
ness 0(x-|) and using the fact that 0(x)/il(x) = 0* = 1 .2386, we have 


1/2 


0 (X) 

0 (x-i ) 


(0/5,)' 


X - x-| 


Ui 0 (XT )1 0 (XT ) 


V 


+ 1 


(Laminar flow) (17) 


As mentioned before, we will assume that the functional form of the mean 
velocity profile in a self-similar turbulent shear layer is the same as in 
laminar flow. Thus, we continue to use equations (12) and (13) for the velocity 
profile, though the definition of the transverse length scale ^ {x) (which 
equals (vx/Ut)^/^ in laminar flow) must be altered. Since the equation 
0=0* £ (x) holds in general, this may serve as the definition of (x) in 
turbulent flow, provided that an adequate empirical formula for the development 
of 0 (x) can be found. 


From the data of Sato (ref. 12), such a formula is not hard to develop. 
We may take, for example, in the transitional and fully turbulent regions of 
the flow 


0 (x) 



where, for (x 


0(x) 


's J 


lam 


Xg)/0g = N^., 



extra 
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tanh 


lx - Xg 
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- Nt 

0s 
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(18) 


(19) 
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and where given by equation (17) with x-j = Xg. A plot 

of equation (19) with the numerical values X-|/9g = 34, X2/9s “ 135, and 

= 56 is shown in figure 6 . Also shown au:e the measured values of 
0(x)/0g - [6(x)/0g]iam three different flows in reference 12. The agree- 

ment is seen to be quite reasonable. 

The condition (s - Xg)/0g = N^. represented a transition criterion and 
appears to hold quite generally for two-dimensional laminar separated flows of 
the type under discussion. For example, figure 12 of Freymuth (ref. 19), which 
shows the growth of the streamwise fluctuation velocity component lu'lmax 
a free shear layer downstream of separation for a variety of separation Reynolds 
numbers, indicates that the range of (x - Xg) 0 g in which the disturbances 
grow exponentially is almost completely independent of Reynolds number for 
61 ^ R 0 ^ ^ 334. This range was found to be 0 ^ (x - Xg)/0g ^ 60. Strong non- 
linearity (indicating the onset of transition) appears for (x - Xg)/ 0 g ^ 60, 
which agrees well with the value N^- = 56 used here. A similar independence 
of the onset of transition with Reynolds number can be found in separation 
bubbles. For example. Caster's results (ref. 14, fig. 11) of intermittency 
measurements in seven short separation bubbles in the range 136 ^ R 0 ^ i 394 

show that the onset of transition occurs for 50 ^ (x - Xg)/0g S 80 for all but 
one of the bubbles. Unlike the case of fully attached flew, it therefore 
appears that transition in fully separated flow is relatively easy to predict 
in terms of the critical value of a single parameter, at least in the case of 
ncxninally two-dimensional incompressible flow in the Reynolds number range 
under discussion under conditions of small background turbulence and acoustic 
noise. 


VELCX:iTy DISTRIBUTION ALOIG A STREAMLINE 

pro- 
0 , 


(20) 


mass 

flow across the entire lower half of the shear layer in figure 5, is a function 
of X and may be used to define the (constant) stream function along the 
streamline that just enters the underside of the shear layer at any station x. 
In our model of the separation bubble, we will suppose that the separation 
point X = Xg corresponds to the point where the separation streamline enters 
the underside of the separated shear layer in a bubble in just this way. That 
is, we will take 


4;s = Ug £(Xg) f(-~) = -Ug0g (21) 


The stream function can be obtained from our approximate velocity 
file (13) by integration. Applying the boundary condition = 0 on n = 
we have, from equations (7) and (13), 


Ui £(x) 3 / 

4; = Ui S,(x) f(n) f'(0)<rtanh 

b 2 ] 


1 


bn + tanh ^ ( — 

VA 


Now the value of the stream function at ri = -<», which represents the 
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as the value of the stream function (relative to the shear layer) that repre- 
sents the separation streamline in a bubble. 


The distribution of velocity along the streamline i]; = ipg can now be found 
easily. From equations (20) and (21) (with Ui = Ug) , we have 


tanh 


1 


bn + tanh”^ 


®s 2 b /e\ 1 

9 3 f(0)\A/ 


( 22 ) 


It follows from equation (13) and the identity sech^ (x) = 1 - tanh^ (x) that 
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f (n) 




= f ( 0)-<1 
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-2b0 


* 0 


s 1 
+ — 


3 f’(0) 0 


(23) 


Figure 7 is a plot of equation (23) for various separation Reynolds numbers, 
Equations (17) through (19) have been used to evaluate the quantity O/Og. 


STAGNATION PRESSURE ALONG ijj = ijjg 

Equation (23) , which was calculated from the known behavior of a free shear 
layer in zero pressure gradient, provides a measure of the stagnation pressure 
distribution along the streamline \1^ = \|^g, since in zero pressure gradient, the 
stagnation pressure H(x,4'), defined by 


H(x,\[i) 



(24) 


is simply equal to the dynamic pressure (p/2)u^. 

We recall that in a general shear flow, the stagnation pressure varies 
along a mean streamline according to the approximate equation 

8h 3t 


where T is the combined viscous and turbulent shear stress 


8u 

T = )J — - pu'v' (26) 

3n 
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and where s and n are orthogonal curvilinear coordinates in the plane of 
the flow oriented in the streamwise and cross-stream directions, respectively. 
Equation (23) thus represents the cumulative effect of such cross-stream stress 
gradients in the downstream direction. 

The use of equation (23) to model the actual stagnation pressure distribu- 
tion along \1 j = \}jg in a separation bubble involves the assumption that the 
modification of the shear stress distribution due to the abrupt rise in static 
pressure at the rear of the bubble is small, or, equivalently, that the turbu- 
lent shear stresses are dominated by "memory" effects. Such an approximation 
of memory-dominated shear stresses is familiar in many analyses of turbulent 
flews following entry into a region of severe adverse pressure gradient (see 
especially refs. 20 and 18). The validity of this approximation when applied 
to the bubble problem has yet to be established, however. 

With due regard for this uncertainty, we may take the nondimens ional form 
of equation (24) as our working equation for the distribution of stagnation 
pressure along the streamline ip = \|jg: 


layer 


where the last term on the right-hand side is given by equation (23) . 



-f 

s /plane 


(27) 


CASTER'S BURSTING CRITERION 


M. Caster (ref. 14) 

parameters Ra and the 
s 


has obtained an empirical correlation in terms of the 
pressure gradient parameter P, defined by 


_ Uj. - Us 0g (Uj. - Us)/Ug UgOg 

P = (28) 

^r — ^s ^ (^r ~ ^s^/^s ^ 


that distinguishes between flew conditions permitting a short separation bubble 
and the conditions when a short bubble is not possible (i.e., a "bursting" con- 
dition). Figure 8 shows Caster's so-called "bursting line," which we will 
attempt to derive from equation (27). 

If we define a pressure coefficient c^ relative to the free stream con- 
ditions at the separation point in a bubble, then 
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In terms 


where Cp = ^ “ (U/U^) ^ is the conventional pressure coefficient, 
of c^, equation (27) becomes, on \|^ = 'iJg, 



layer 


At the reattachment point, the streamwise velocity u on 4^ = il;g equals zero, 
so that 


/ Reattachment'' 
\ condition j 




plane 

layer 


(31) 


where the pressure recovery coefficient a is the same as in equation ( 1 ). 

Equation (31) together with equation (28), which at reattachment can be 
written as 


P 


- cgr - 1 

(x Xg ) /9 g 


(32) 


define the largest possible P that can be achieved by a short bubble at a 
given K0g* That is, if Cpj- frcxn equation (31) is substituted into equa- 
tion (32) and the resulting P tabulated as a function of (x - Xg)/0g for a 
particular then an extremum of P is found. For the Reynolds numbers 

R 0 g =100, 200, 300, and 400, we find that the extreme values of P are 

-0.0996, -0.1868, -0.2732, and -0.3593, occurring at (x - Xg)/9g =103, 111, 
115, and 117, respectively. A faired curve through these points is shown in 
figure 8 and exhibits suprisingly good qualitative agreement with Caster's 
empirical curve. 

The agreement between theory and experiment is particularly pleasing in 
view of the fact that all of the empirical input to the present calculation 
method was obtained from measurements of transition in free shear layers in 
zero pressure gradient at a remote distance from any wall. Our calculations 
therefore lend support to the view that the processes of bubble bursting and 
shear layer reattachment are linked to the processes of laminar-turbulent 
transition and shear layer growth in a fundamental way. 
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CDMPARISCN WITH EXPERIMENTAL BUBBLE LENGTHS 

As a caution against overly optimistic conclusions drawn from the present 
analysis, we present in figure 9 the results of an attempt to verify the 
reattachment condition (31) directly. In the figure, we have plotted the mea- 
sured pressure recovery coefficient o and length (Xj. - Xg)/6g for a number 
of separation bubbles from various sources. On the same scale, we have drawn a 
copy of figure 7, which represents the nondimens ional stagnation pressure along 
the separation streamline for a plane layer. If equation (31) held exactly, 
then all of the experimental points would lie on the dashed curves. 

In fact, the experimental points are scattered quite widely about the 
analytic curves. More serious than the scatter in o, however, is the discrep- 
ancy between the measured length of the bubbles, which are typically in the 
range 150 ^ (x^ - Xg)/0g ^ 300, and the values calculated from our analytic 
derivation of Caster's bursting line. We have indicated the latter points in 
figure 9 by "plus" symbols, which lie, to within graphical accuracy, on the 
family of straight lines drawn through the origin and the points of tangency 

2 

with the various curves of (u/U){j)-\|j^ for the plane layer. As can be seen, 

the calculated bubble lengths are typically a factor two or three too small. 

This error is significant, and indicates that the present theory is not suf- 
ficiently refined for practical use in the quantitative prediction of bubble 
lengths . 

It is not clear which of the many effects neglected in the theory are 
responsible for this discrepancy. We may speculate, however, that the neglect 
of changes in the shear stress gradient Sx/Sn on the separation streamline 
during reattachment is important. Seme such change must occur, since in the 
real flow the shear stress on \|) = vanishes at reattachment, which places 

the reattachment point belcw the point of maximum T in the layer. The crude 
theory developed above assumes that as the separation streamline gets swallowed 
up into the underside of the shear layer, the shear stresses on the streamline 
tend asymptotically toward the value of the center of the layer (which is the 
maximum value) . It follows that 3x/3n is underestimated by the crude theory. 
This implies (from equation (25)) an unrealistically low prediction of the stag- 
nation pressure H. It is not obvious how the present theory could be refined 
to take such effects into account. 


CONCLUSIONS 

From our investigation of the physics of separation bubbles, we may state 
the following tentative conclusions; 

1 . The form of the static pressure distribution about a short separation 
bubble can be explained qualitatively in terms of boundary layer concepts. In 
particular, the sharp corner in the pressure distribution at the rear of a short 
bubble is more properly interpreted as the point where the lower edge of the 
reattaching shear layer grazes the surface than as the point of transition frcrni 
laminar to turbulent flow, as has traditionally been assumed. 
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2. The observation that the static pressure distribution about the body 
on which a short separation bubble occurs very nearly coincides with the 
inviscid pressure distribution at all points other than those between separa- 
ticxi and reattachment leads to a plausible hypothesis regarding the factors 
controlling bubble bursting and a simple but crude calculation method for esti- 
mating the length of a short bubble when it is possible for one to exist. Spe- 
cifically, the reattachment point coincides with the earliest downstream 

statical X at which the stagnation pressure on the separation streamline is 
large enough to support an abrupt rise to the local inviscid static pressure 
at X. 


3. The attempt to approximate the rise in stagnation pressure on the sepa- 
ration streamline in a bubble by the corresponding rise along a streamline 
entrained into the underside of a laminar mixing layer undergoing transition 
leads to a reasonably accurate prediction of Caster's bursting line, but signif- 
icantly underestimates the length of short bubbles in general. This inaccuracy 
is probably due to the failure of the calculation method to take account of 
changes in the viscous and turbulent shear stress profiles occurring prior to 
reattachment that produce (according to equation (25)) an increase in the actual 
stagnation pressure along the separation streamline beyond what exists in a free 
shear layer. 

4. The present results are in every way consistent with and support the 
view of Cebeci and Bradshaw and others that reattachment is a direct consequence 
of the increased entrainment of fluid into the underside of the growing sepa- 
rated shear layer downstream of transition, the entrainment causing the shear 
layer to suck itself back down onto the surface. 
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APPENDIX 


We rewrite equation (15) in the equivalent form 


rtao 

X E’- 


f) - (1 - f) 




(f ')2 dn 


(A1) 


Substituting equations (12) and (13) into equation (A1 ) gives 

|[l - f (0)]li - [l - f (0)J^I^1 = I3 


(A2) 


where 


■■ ■ r& r 


0 \Jt[ '^ari /2 


e dx) a dn = -;= 


f. 


(A3) 


■ X fs X„„ "■') " 


^ yfir ^^ 311/2 


]fii 


(A4) 


and 


fO (3 

00 )2 


SGCh^ 


bn + tanh ^ j — 

\/3> 


b dn = 


4J3 - 9 


(A5) 


(assume b < 0 throughout) 

The right-hand side of equation (A2) can be rewritten by means of equation (14) 
as 


I 3 f ' (0) 2 ^ _f, ( 0 ) 3 2 /rr 1 

b ” 1 - f ( 0 ) ^ a 


(A 6 ) 


Substituting the identity 


-f • (0) 3 = [1 - f ' (0)] ^ - 3 [1 - f (0)J ^ + 3(^1 - f • (0)] - 1 
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into equation (A 6 ) and combining with equation (A2) gives the cubic equation 

[l - f-( 0 )J 3^-12 - I 3 + D - f'(0)]2^Il + 3 I 3 


/ 2|/i\ 2 i/tt 

[l - f • (0)] \-2l2 — j + I 3 — = 0 


(A7) 


which may be solved numerically by Newton's method to give the solution 
1 - f'(0) = 0.41077, or f ' (0) = 0.58923. 

The constants a and b may be calculated from the normalization condi- 
tion (16) in the following way. The first integral on the right-hand side of 
equation (16) may be written 


fO r f ( 0 ) f'( 0)2 

Lf - <3n = I 4 

«^_oo b b 


-f'( 0 ) 2 |/i 


|_1 - f ( 0 ) 


^^'(0) I 4 - [f'(0)J^I^ (A 8 ) 


where we have used equation (14), and where 


= r° - sech2 
^-co 2 


bri + tanh“^ / — ' 

\fv 


b dn = 


J 3 - 3 


(A9) 


The second integral on the right-hand side of equation (16) is just the left- 
hand side of equation (A2) . According to equation (16), the sum of these two 
integrals equals -f(-°°) = 1 .2386, which immediately gives a = 0.88544. Equa- 
tion (14) then gives b = -0.30160. 
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ACTUAL PRESSURE DIST. 


INVISCID PRESSURE DIST. 




Figure 1.- Characteristic pressure distribution beneath a 
separation bubble of the "short" type. 


EDGE OF SHEAR LAYER 
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Figure 2.- Two-dimensional steady separation from a 
solid wall showing entrainment of fluid into the 
underside of the shear layer due to spreading in 
the downstream direction. 
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Figure 3.- Recirculating flew pattern within a bubble 
showing splitting of the shear layer at reattach- 
ment to supply the entrainment (vertical scale 
greatly exaggerated) . 


Figure 4.- Pressure distribution inferred fron figure 3. 
Small entrainment velocities beneath shear layer imply 
small pressure gradients between S and K. Steep 
gradient between K and R due to splitting of shear 
layer at reattachment and gradient of stagnation pres- 
sure across lower half of shear layer . 





DIVIDING / ' 

STREAMLINE 


Figure 5.- 1 


FLOW 

RELATIVE 
GROWTH IN 
MOMENTUM 
THICKNESS 


1 5 ^ 


1 0 


0 




0 


T 

TRANS 



SQUARED 
VELOCITY 
ON y= Vg 

IN ZERO 

PRESSURE 

GRADIENT 


Figure 


STRENGTH 

ADVERSE 

PRESSURE 

GRADIENT 



7.- Plot of equation (23) using equations (17), (18), and (19). 
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Figure 8.- Gaster's bursting line derived from present theory 
which is based entirely on data from free shear layers. 




Figure 9.- Attempt to verify reattachment condition (31). 
Measured lengths and pressure recovery coefficients 
frcm various sources: q. Caster (ref. 14) "series I"; 

A, Caster "series II"; Or Van Ingen (ref. 11); V, Tani 
(ref. 2, fig. 17, based on data frcxn Cault, ref. 10), 

If equation (31) held exactly, then all measured points 
should lie on dashed curves. Note: each set of data 

represents a series of runs at various tunnel speeds, 
which, to avoid cluttering the figure, have not been 
shown. 



